How stable are transport model results to changes of resonance parameters? 

A UrQMD model study 



Jochen GerhardQ Bj0rn Bauchle, Volker Lindenstruth, and Marcus Bleicher 
Frankfurt Institute for Advanced Studies, Ruth-Moufang-Strafie 1, 60438 Frankfurt am Main 
Institut fur Theoretische Physik, Johann Wolfgang Goethe-Universitdt, 
Max-von-Laue-Strafle 1, 60438 Frankfurt am Main and 
Institut fur Informatik, Johann Wolfgang Goethe-Universitdt, 
Robert- M ay er-Strafie 11-15, 60054 Frankfurt am Main 

The Ultrarelativistic Quantum Molecular Dynamics [UrQMD] model is widely used to simulate 
heavy ion collisions in broad energy ranges. It consists of various components to implement the 
different physical processes underlying the transport approach. A major building block are the 
shared tables of constants, implementing the baryon masses and widths. Unfortunately, many of 
these input parameters are not well known experimentally. In view of the upcoming physics program 
at FAIR, it is therefore of fundamental interest to explore the stability of the model results when 
these parameters are varied. We perform a systematic variation of particle masses and widths 
within the limits proposed by the particle data group (or up to 10%). We find that the model 
results do only weakly depend on the variation of these input parameters. Thus, we conclude that 
the present implementation is stable with respect to the modification of not yet well specified particle 
parameters. 



I. INTRODUCTION 

One of the major themes of todays high energy 
physics is the exploration of Quantum Chromo-Dynamics 
(QCD). QCD predicts that at sufficiently high tempera- 
tures and densities, nuclear matter could exhibit a phase 
transition into a new state of matter the Quark Gluon- 
Plasma (QGP). In this state the usual color confinement 
is relaxed and the constituents of the matter, namely 
quarks and gluons, are allowed to move over distances 
larger than the scale of a single nucleon. Indeed, experi- 
ments at the CERN-SPS, the BNL-RHIC and the LHC at 
CERN have collected an impressively large body of data 
that is consistent with the interpretation that a QGP 
was formed for a short period of time. While the current 
LHC program is running at the highest available ener- 
gies and therefore provides insights into the properties 
of the QGP at high temperatures and very low baryon 
densities, the RHIC program is now focused on the ex- 
ploration of the QGP phase transition with a low energy 
scan program. Here one hopes to find the existence and 
location of the critical end point of the first order transi- 
tion line that is expected at high baryon densities. From 
the year 2018 on, this scan for the onset of deconfinement 
and the search for the critical end point will also be a top 
priority of the FAIR facility currently under construction 
in Germany. 

Apart from the experimental difficulties, a major ob- 
stacle to pin down the properties of strongly interact- 
ing matter is the unambiguous interpretation of the ex- 
perimental results. Unfortunately, first principle lattice 
QCD calculations are currently only feasible in thermal 
equilibrium and for very moderate T//ib values (T be- 
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ing the temperature and \xb being the baryo-chemical 
potential). Therefore, transport approaches like the 
Ultra-relativistic Quantum Molecular Dynamics model 
(UrQMD), the Parton-Hadron-String-Dynamics (PHSD) 
PQ, the Multi-Phase- Transport model (AMPT) and 
many other dynamical models are employed to link the 
final state observables to the physics properties of the hot 
and dense stage of the reaction. All these models have in 
common that they rely as input on measured quantities 
like the hadron masses, the hadron decay widths, individ- 
ual branching ratios and cross sections. Unfortunately, 
these quantities are very often not exactly known, as one 
can see from an inspection of the Particle Data Group 
(PDG) tables. 

In this paper we explore systematically how a varia- 
tion of some of these parameters influences the results 
obtained in transport simulation. As an example, we 
employ the UrQMD model[MS]- In the long run, these 
investigations will allow to obtain a systematic error of 
the simulations, which is needed to quantify the quality 
of the model results. While we restrict ourselves to the 
UrQMD model in this study, the results should (qualita- 
tively) also be transferable to other transport simulations 
based on similar physics assumptions as the ones men- 
tioned above. 



II. COMPUTATIONAL SET-UP 

The idea is to replace the hard-coded hadron masses 
and widths with automatically generated tables with var- 
ied parameters. Then a sufficient number of simulations 
is performed and evaluated. To this aim, we have de- 
signed Python modules to read the up-to-date data from 
PDG web page [6]. The usage of Python as glue code for 
high performance computations allows for very short de- 
velopment cycles by combining rather slow, but powerful 
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FIG. 1. Top: Total pion-nucleon cross section as a function of 
^/s for a systematic variation of the nucleon resonance masses. 
The full line depicts the PDG averages, the dotted line shows 
a variation of the PDG parameters by +10%, the dashed line a 
variation by -10%. Bottom: Total pion-nucleon cross section 
as a function of -^/s for a systematic variation of the nucleon 
resonance widths. The full line depicts the PDG averages, 
the dotted line shows a variation of the PDG parameters by 
+10%, the dashed line a variation by -10%. 



scripting parts of the program with very fast compiled 
parts in another language. This possibility is often used 
in graphics card computations [7 . We use Python to au- 
tomatically rewrite the Fortran source code of UrQMD 
according to recalculated variations of the PDG-Data. 
This form of meta programing is often preferable to us- 
ing configuration files. With the parameters being hard- 
coded in the Fortran source they are known at compile 
time, which enables the compiler to do more optimiza- 
tions. We employ the Frankfurt LOEWE-CSC[8] to carry 
out a systematical parameter scan and to check of the 
stability of UrQMD. After the computation the data files 
from different UrQMD runs are parsed and compressed 
to statistics files, which can easily be interpreted on local 
systems. 

We explore the dependence of the UrQMD results on 
the particle data, within the estimated errors provided 
by the PDG. In addition we vary the parameters within 
an overall range of ±10% of mass and width. Separate 



scans for variations of mass and width of each baryon 
family are performed. Typically up to lO'OOO events are 
simulated to stay clear of statistical errors. 



III. CROSS SECTIONS 

Let us start by investigating the total pion-nucleon 
cross section as a function of energy. This cross sec- 
tion is of special importance for the dynamics of nuclear 
matter at intermediate energies. It also serves as direct 
benchmark to adjust the parameter sets since it is well 
measured experimentally. The total cross section is given 
by 

<TnI= ^2 tiN,m N ,j 7T ,m 7r \\JB.,M R )x 

%Sr + 1 7T rij^Ar^rtot /-.n 

X (2S N + l)(25 ff + 1) p 2 CMS ( Mfl _ ^)2 + Ik ' 

with the total and partial decay widths Ttot and Tn^N n . 
Thus, the cross-section depends on the widths and masses 
of all nucleon- and delta-resonances N* and A^*). Fig- 
ure [T] (top) depicts the total pion-nucleon cross section 
cr£j? as a function of the center of mass energy i/s for 
a systematic variation of the nucleon resonance masses. 
The full line depicts the PDG averages, the dotted line 
shows a variation of the PDG parameters by +10%, the 
dashed line a variation by -10%. One clearly observes 
that the 7tN cross section varies strongly if the resonance 
masses are changed. In turn, however, this allows to pin 
down the resonance masses rather precisely. In contrast 
a change of the resonance widths leaves the cross sec- 
tion unaltered. Figure [I] (bottom) shows the total pion- 
nucleon cross section as a function of ^fs for a systematic 
variation of the nucleon resonance widths. The full line 
depicts the PDG averages, the dotted line shows a vari- 
ation of the PDG parameters by +10%, the dashed line 
a variation by -10%. 

Figure [2] (top) depicts the total pion-nucleon cross sec- 
tion as a function of y/s for a systematic variation of the 
A masses. The full line depicts the PDG averages, the 
dotted line shows a variation of the PDG parameters by 
+10%, the dashed line a variation by -10%. Figure [2] 
(bottom) shows the total pion-nucleon cross section as a 
function of y/s for a systematic variation of the A widths. 
The full line depicts the PDG averages, the dotted line 
shows a variation of the PDG parameters by +10%, the 
dashed line a variation by -10%. One clearly observes 
that a variation of the nucleon masses has a drastic ef- 
fect on the pion-nucleon cross sections especially in the 
A(1232) region. In comparison to the available experi- 
mental data, strong constraints on the model parameters 
can be obtained. In fact, the employed parameters are 
based on the PDG data and re-adjusted within the limits 
of the PDG ranges. 
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FIG. 2. Top: Total pion-nucleon cross section as a function of 
yfs for a systematic variation of the A masses. The full line 
depicts the PDG averages, the dotted line shows a variation of 
the PDG parameters by +10%, the dashed line a variation by 
-10%. Bottom: Total pion-nucleon cross section as a function 
of \fs for a systematic variation of the A widths. The full line 
depicts the PDG averages, the dotted line shows a variation 
of the PDG parameters by +10%, the dashed line a variation 
by -10%. 



IV. PION PRODUCTION 
A. Pion yields 

Let us next turn to the investigation of full Pb+Pb col- 
lisions and focus on the FAIR energy range of 2 AGeV 
and 30 AGeV. Here we investigate the total pion yield 
for a systematic variation of all nucleon resonance masses 
mjv* by up to 10%. We show the deviation of the pion 
yield compared to a UrQMD calculation with the mean 
values of the PDG-Data files. Figure [3] (top) shows 
the relative pion yield in Pb+Pb collisions at 2 AGeV 
beam energy for a systematic variation of the masses 
and widths of the nucleon resonances by ±10%. Fig- 
ure [3] (bottom) shows the pion yield in Pb+Pb collisions 
at 30 AGeV beam energy for a systematic variation of the 
masses and widths of the nucleon resonances by ±10%. 

Even at the lowest energy, which is strongly dominated 
by resonance dynamics, the model results do at worst 



FIG. 3. Top: Relative pion yield in Pb+Pb collisions at 
2 AGeV beam energy. For a systematic variation of the 
masses and widths of the nucleon resonances by ±10%. Bot- 
tom: Relative pion yield in Pb+Pb collisions at 30 AGeV 
beam energy. For a systematic variation of the masses and 
widths of the nucleon resonances by ±10%. 



vary linearly with the variation of the model parame- 
ters. At 30 AGeV, the model results are stable against a 
variation of the resonance parameters. The variation of 
the particle widths has no significant effect on the model 
results. 

Figure [^investigates the pion production as a function 
of varying masses of the delta-resonances jtia and their 
widths Ta ■ The masses of all delta resonances have been 
scaled with the same factor. Here we limit the varia- 
tion to -8% — 1-10%, because the code becomes unsta- 
ble for too low masses. Again a variation of the width 
leaves the results unchanged. The variation of the A^*) 
masses, however results in a strong variation of the pion 
yield at 2 AGeV. This effect is mainly attributed to the 
A(1232) resonance that is pushed towards the kinematic 
limit (mA(i232) ~~ ^ m p ± fn^)- At 30 AGeV, the variance 
of the yield stays generally moderate. However, the pion 
yield shows a pronounced step if the masses are shifted 
by ~ 1.8% upwards. While the magnitude of the effect 
is small it indicates that complex simulation models may 
exhibit discontinuous behaviors. Let us now have a closer 
look into the origin of the step. 
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FIG. 4. Top: Relative pion yield in Pb+Pb collisions at 
2 AGeV beam energy. For a systematic variation of the 
masses and widths of the A resonances by ±10%. Bottom: 
Relative pion yield in Pb+Pb collisions at 30 AGeV beam 
energy. For a systematic variation of the masses and widths 
of the A resonances by ±10%. 



For further analyses, we group different production 
processes into 5 classes, discriminating the decay of delta- 
resonances (A), the decay of nucleon-resonances (N*), 
the decay of strange baryons (accounting for all un- 
stable baryons not included in the former two classes) 
(B s ), the decay of meson resonances (m) and scatterings 
(XY -± it + R). In Figure [| (top) one observes that at 
E\a,b = 2 AGeV, the number of pions from scatterings 
stays constant as a function of delta mass, while the pro- 
duction of pions from delta decays rises linearly. At very 
low delta masses (less than -6% of the standard value), 
the formation of deltas absorbs on average more pions 
than are being produced by the decays thereof, while at 
higher delta masses, the opposite is true. The increase of 
pion production from deltas is counteracted by a decrease 
of pion production from nucleon resonances N* , which 
start to be net-absorbing above ±6% of the standard 
(delta-)masses. This can be explained in a picture of de- 
tailed balance: When the delta resonances produce more 
pions, the equilibrium value of pion- and N* -multiplicity 
is shifted towards the N* . Thus, the iV*-phase space 
is populated more quickly than it is depleted. The same 



FIG. 5. Top: Pion yield in Pb±Pb collisions at 2 AGeV 
beam energy itemizing different production processes. Bot- 
tom: Pion production from various A resonances in Pb±Pb 
collisions at 30 AGeV. 



effect, though much weaker and not turning around com- 
pletely, can be seen in the decrease of the number of pi- 
ons from mesonic decays. In total, the rise of pions from 
deltas counteracts the fall of pions from N* , thus leading 
to a weak overall rise of the pion production. 

At 30 AGeV we focus now on the step like behavior 
at a A-mass shift of ~ ±1.8%. Figure [5] (bottom) shows 
the contributions of different A-resonances to the final 
number of pions for varying delta masses between and 
±10% at high impact energy E\ a b — 30 AGeV. I an analy- 
sis simlar to the one from Figure[5](top), we trace the step 
to the delta contribution. The step we discovered earlier 
consists of an increased pion production (less absorption) 
from the Aig5o-resonance, which rises from -20 to 0, and 
a corresponding decreased production (increased net ab- 
sorption) from the Aig 2 o-resonance, which drops from 
to -16. The y^-distribution of the underlying ttN colli- 
sion remains essentially unchanged as a function of itia- 
Therefore, at the point of the discontinuities, the A 192 o 
takes the role that Aig5o had at lower masses. Since the 
branching ratios Ai 92 o -> 7rAi 2 32 and A1950 -> 7rAi 232 
differ by a factor of 2 (40% vs. 20%), the number of A i232 
and thus the number of pions are changed over a small 
mass intervall. We find a simliar behavior at ~ ±5%, 
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FIG. 6. Top: Deviation of pion transverse momentum spec- 
tra (pt) in Pb+Pb collisions at 2 AGeV for variations of the 
nucleon resonance masses. Bottom: Deviation of pion trans- 
verse momentum spectra (pt) in Pb+Pb collisions at 30 AGeV 
for variations of the nucleon resonance masses. 



FIG. 7. Top: Deviation of pion transverse momentum spec- 
tra (pt) in Pb+Pb collisions at 2 AGeV for variations of the 
A resonance masses. Bottom: Deviation of pion transverse 
momentum spectra (pt) in Pb+Pb collisions at 30 AGeV for 
variations of the A resonance masses. 



where the pion production from Aigso decreases, while 
the production from Aigoo increases. Superimposed is 
an approximately linear rise of pion production from the 
lowest delta resonance A 12 32- 



B. pT-spectra of tt + 

Finally, we discuss the transverse momentum 
distributions [9] ■ Here we investigate the pions transverse 
momentum distributions in Pb+Pb collisions at 2 and 
30 AGeV. Figure [6] (top) shows the deviation of the pion 
transverse momentum spectra (pt) in Pb+Pb collisions at 
2 AGeV for variations of the nucleon resonance masses. 

At 2 AGeV, a decrease of the nucleon resonance masses 
shifts the pions to lower p t , while an increase of the 
masses shifts it to higher transverse momenta. However, 
variations of the yields at higher p t maybe up to ±20% 
for a variation of ±5%. 

Figure [6] (bottom) displays the deviation of the pion 
transverse momentum spectra (pt) in Pb+Pb collisions 
at 30 AGeV for variations of the nucleon resonance 
masses. Figure [7] (top) shows the deviation of the pion 



transverse momentum spectra (pt) in Pb+Pb collisions 
at 2 AGeV for variations of the A resonance masses. 
Again, we do not observe and effect on the calculations 
at E\ a b = 30 AGeV is within the statistical fluctuations. 
Figure [7] (bottom) displays the deviation of the pion 
transverse momentum spectra (pt) in Pb+Pb collisions 
at 30 AGeV for variations of the A resonance masses. 

The effect of variations in delta resonance masses is 
strongly non- linear. Both at high beam energies and at 
low beam energies, we can distinguish three transverse 
momentum regions. Pions from intermediate transverse 
momentum 0.2 < p t < 0.6 GeV are being shifted to low 
transverse momentum p t < 0.2 GeV, if the masses are 
decreased and vice versa, if the masses are increased. 
This is expected, since the available kinetic energy in 
a delta decay decreases with decreasing delta mass. At 
higher transverse momenta, lower masses lead to higher 
pion yields, while higher masses lead to lower pion yields, 
which is a reversal from the behavior observed from vary- 
ing the nucleon resonance masses. Furthermore, we ob- 
serve the effects to be a lot stronger in low-energy colli- 
sions. Also the variation of the A masses by ±5% results 
in modifications of the pion yield by ±20% in given p t 



G 



regions. 

V. SUMMARY 

In light of the upcoming high precision experiments at 
FAIR, it is highly desirable to obtain better estimates 
on the systematic errors of transport simulations. We 
addressed this question by using the UrQMD transport 
approach in nucleus-nucleus reactions in the FAIR energy 
regime from 2 to 30 AGeV. We have analysed elementary 
cross sections in pion-nucleon reactions, as well as lead- 
lead collisions for various sets of input parameter varia- 
tions of the hadron masses and widths. Although the an- 
alyzed quantities show globally only a weak dependence, 
discontinuities like in Figure [t] (bottom) may occur and 
influence predictions made by the applied models. The 
dependence is strongest in low-energy collisions, where 
the collision dynamics is dominated by resonance pro- 



duction and decay. Here, one may encounter systematic 
errors on the order of ±20%. At higher energies, the sys- 
tematic errors are much smaller. One should note, that 
the present study explored a worst case scenario where 
all parameters were shifted simultaneously in one direc- 
tion. The error on the masses (and widths) are however 
uncorrelated and should therefore induce smaller system- 
atic bias into the simulations as compared to this study. 
Therefore, we conclude that the predictive and analysis 
power of the present approach is better than a systematic 
error of 20 %. 
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